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1  SUMMARY  OF  ACCOMPLISHMENTS 


1.1  Characteristics  of  Program 

First  ever  demonstration  of  military  critical  capabilities; 

•  High  dynamic  range  transient  simulation 

•  First  full-wave  interconnect  Modeling  Tool 

•  Only  way  to  model  critical  technologies 

•  Molecular  electronics 

•  Tuimeling  in  nanoelectronics 

•  Radar  systems  (transient) 

•  Driver  of  commercial  developments 

•  Open  Source 

1.2  Software  tools 

We  have  developed  five  software  packages  responding  to  feedback  from  users. 
These  packages  are: 


•  fREEDA 

o  High  dynamic  range  multi  physics  simulator 
o  Easy  development  of  advanced  device  models 
o  Many  device  models  never  before  implemented 

•  S2IBIS3  (Spice  To  IBIS) 

o  Automated  digital  macromodeling  tool 

•  EMPDK 
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o  EM- Aware  physical  design  kit  tool  (JAVA) 

•  UIUC2D 

o  EM  modeling  tool  for  2D  geometries 

•  ICWAVE 

o  On  chip,  comprehensive,  3D  EM  modeling 
o  Technology  used  in  part  in  NeoRF  (NeoLinear) 

•  PRIME 

o  EM  Reduced  Order  Model  interconnect  modeler 
o  Directly  interfaces  with  fREEDA 
o  Synthesizes  R,  L,  C,  K  models  for  Spice. 

All  packages  can  be  accessed  through  http://www.freeda.org 
PRIME  and  UIUC2D  http :  //alpha  1 .  ece .  uiuc .  edu/download 
For  PRIME:  username:  prime  pwd:  fitting 
For  UIUC2D:  username:  uiuc2d  pwd:  rlcgsyn 
For  ICWAVE  Contact  Andreas  Cangellaris 

We  will  continue  to  maintain  these  products. 

1.3  Modeling  of  Mixed  Signal  X-Band  Radar 

Traditional  microwave  subsystem  design  is  based  on  harmonic  balance  (frequency-domain) 
simulations.  However  real-world  effects  due  to  transients,  e.g.  temperature  during  pulsed 
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Transmission  are  not  captured.  Many  effects  are  signal  dependent  and  the  chirp-like  radar  signal 
cannot  be  represented  in  the  frequency  domain.  Only  a  transient  simulation  will  capture  all 
effects  but  the  problem  here  is  that  transient  simulation  has  typically  low  dynamic  range  and  the 
undesirable  effects  such  as  spectral  regrowth  cannot  be  seen.  Many  of  the  developments  required 
to  develop  a  transient  simulator  with  all  of  the  attributed  required  for  circuit-level  radar 
simulation  are  described  in  subsequent  sections. 

The  canonical  radar  transmit  system  shown  in  Figure  1  was  modeled.  The  mixer  was  considered 
to  be  ideal  but  a  high  order  bandpass  filter,  a  circuit-level  model  of  an  X-band  MMIC,  and  a 
behavioral  model  of  a  Traveling  Wave  Tube  Amplifier  (TWTA)  were  modeled.  The  TWTA 
model  requires  long  time  delays  that  have  never  been  modeled  in  a  transient  simulator 
previously. 

Figure  2  shows  the  transient  waveform  at  the  output  of  the  TWTA.  The  initial  spike  results  from 
turning  the  system  on.  This  is  clearly  an  effect  that  must  be  suppressed  and  its  level  could  not  be 
determined  in  a  traditional  harmonic  balance  simulation.  The  spectrum  at  the  output  is  shown  in 
Figure  3.  It  was  not  possible  to  use  an  FFT  to  arrive  at  this  spectrum  with  this  fidelity.  Instead 
first  principles  Fourier  transformation  was  used.  We  believe  that  these  results  are  the  first 
circuit-level  transient  simulation  of  a  radar.  We  are  currently  constructing  this  radar  but 
comparisons  between  simulations  and  measurements  are  not  available  for  this  report.  Note  that 
actual  bandpass  filter,  driver  amplifier  and  TWTA  were  modeled  and  models  of  the  individual 
subsystems  have  been  experimentally  verified. 
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Figure  2.  Transient  waveform  at  the  output  of  the  TWTA. 
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radar  7. f  ft  .out 


Figure  3.  Spectrum  at  the  output  of  the  TWTA. 


1.4  Multi  Physics  Capability 


fREEDA  is  the  only  simulator  with  advanced  device  models  capturing  full  physics.  Examples  are 
models  of  tunneling  in: 


•  Tunnel  Diode 

•  MOSCAP,  See  Figure  4 

•  TFT 

•  Molecular  Diode,  See  Figure  5 
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The  tunnel  process  cannot  be  expressed  as  a  current-charge-voltage  expression.  We  must  use 
state  variables  to  implement  tuimeling  correctly.  In  addition  fREEDA  is  the  only  simulator  that 
can  handle  long  time  delays  of  more  than  2  or  3  time  steps.  Figure  6  shows  the  input  and  output 
waveforms  of  the  X-band  TWTA. 


Figure  4.  Molecular  Electronics 


Figure  5.  MOSCAP 
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Figure  6.  Input  and  output  waveforms  of  the  TWTA  demonstrating  the  ability  to  model  long  time 

delays. 


1.6  Benchmark  Data 

NeoCAD  benchmark  data  (microwave  measurements)  developed  in  first  year  of  NeoCAD.  Used 
by  EM  modeling  groups  to  validate  data  and  used  in  several  NeoCAD  publications.  Only 
measured  on-chip  benchmark  characterizations  openly  available.  The  three  main  categories  of 
devices  benchmarked  are  sown  in  Figures  7,  8  and  9. 
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Figure  7.  X-Band  MMIC 


(a) 


(b) 


Figure  8.  Silicon  interconnect  characterization  (0.25  Micron). 
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Figure  9.  Coupled  Inductor  (0.25  micron). 


1. 7  NCSU/NeoCAD  Technology  Transfer 

fREEDA  (multiphysics  circuit  simulator) 

•  Open  source,  anyone  can  study  technology. 

•  A  PhD  graduate  of  the  program  moved  to  Synopysis  where  he  led  part  of  the  technical 
development  of  their  latest  simulator  based  in  part  on  fREEDA  Technology. 

•  A  senior  executive  at  Agilent  told  me  that  fREEDA  technology  has  been  implemented  in 
their  ADS  simulator. 

•  A  developer  at  Cadence  Design  Systems  has  indicated  that  fREEDA  has  been  studied 
extensively  in  the  Spectre  group  at  Cadence. 

ICWAVE  (Full  wave  on-chip  EM  simulator) 

•  Vladimir  Okhmatovski,  PostDoc  Univ.  Illinois,  moved  to  NeoLinear  and  worked  on  the 
development  of  NeoRF  partly  incorporating  ICWAVE  Technology 

S2IBIS  (Automated  digital  macromodel  tool) 

•  Extremely  successful,  major  reason  IBIS  macromodels  have  been  successful. 

•  Widely  used:  10  downloads  per  day.  Support  provided  by  EIA  consortium.  Used  by 
virtually  all  IC  vendors.  Impacted  possibly  every  military  electronic  systems 


EM-PDK 


•  EM- Aware  Physical  Design  Kit 

•  Used  by  several  NeoCAD  groups  to  process  interconnect  data. 


The  role  of  a  university  in  computer  aided  engineering  is  to  show  what  is  possible  and  facilitate 
the  transfer  of  tools  and  know-how  to  commercial  CAD  companies.  The  most  successful  method 
of  transfer  is  students  transferring  to  companies. 
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2  Details  on  fREEDA  Development 

2.1  fREEDA  VERSUS  SPICE 

1 .  fREEDA  is  a  multi  physics  simulator 

(a)  It  can  do  things  that  cannot  be  done  any  other  way. 

(b)  High  transient  dynamic  range. 

(c)  Mixed  physics;  Long  time  steps;  Rapid  model  development. 

(d)  Mixed  Electro  /  Thermal  /  EM 

(e)  Many  first  open  source  implementations  of  models. 

(f)  Advanced  concepts,  automatic  differentiation,  general  purpose. 

(g)  Serves  as  a  demonstration  of  what  can  be  done. 

(h)  As  open  source  it  is  being  copied  and  concepts  utilized  by  many  companies. 

2.  HOWEVER,  it  is  not  a  Spice  replacement. 

(a)  E.G.  Spice  and  fREEDA  solve  a  nonlinear  problem  differently. 

(b)  If  the  device  model  is  not  physically  consistent  then  different  solvers  will  give 
different  results.  (This  is  the  case  with  BSIM  models) 

(c)  Works  for  all  other  models,  (e.g.  EKV) 

2.2  Time  Step  Development 

A  major  technology  that  limits  the  dynamic  range  of  a  circuit  simulator  is  the  choice  of  time  step 
and  the  associated  estimate  of  error.  A  transient  simulator  must  estimate  its  error  at  each  time 
step  and  adjust  its  time  point  selection  accordingly.  In  Spice  the  algorithm  is  to  compare  the 
nonlinear  solution  to  a  linear  projection  from  the  solution  at  the  last  time  point.  One  side  effect 
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of  this  is  that  the  time  step  shrinks  as  the  derivative  of  a  waveform  changes  sign,  for  example 
over  the  top  of  a  sinewave.  More  time  points  than  are  needed  lead  to  additional  numerical  noise. 
The  other  fact  is  that  there  is  a  better  way  to  estimate  error.  In  fREEDA  two  different  nonlinear 
solutions  are  performed  and  these  are  compared,  see  Figure  10.  This  is  a  much  better  estimate  of 
error  especially  “over  the  top  of  a  sinewave.”  The  result  is  that  a  better  time  step  is  used  and 
higher  dynamic  range. 


Result  from 
Backward 
Euler 


Result  from 


Error 


Tolerance 


(b) 

Figure  10.  Estimation  of  error  in  fREEDA;  (a)  time  point  selection;  and  (b)  effect  of 

specification  of  tolerance  on  error. 
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2.3  Validation  of  dynamic  range 


Various  definitions  could  be  used  for  the  dynamic  range  of  a  simulator.  The  bets  definition  to  be 
used  should  be  based  on  the  types  of  signals  to  be  captured.  In  the  NeoCAD  project  we 
proposed  the  adoption  of  the  metric  of  the  smallest  detectable  signal  (a  single  frequency  tone)  in 
the  presence  of  a  tone  close  in  frequency.  This  metric  was  adopted  by  others  in  the  NeoCAD 
community.  Data  supporting  our  claim  is  shown  in  Figure  1 1 .  In  Figure  1 1(a)  is  the  layout  view 
of  a  High  Dynamic  Range  PHEMT  2-stage  Low  Noise  Amplifier  MMIC  operating  from  8.5 
GHz  to  14  GHz  with  18dB  Gain.  It  has  a  +6  V  Supply  Voltage,  2  dB  Noise  Figure.  It  is  used  as 
a  pre-driver  amplifier  for  European  phased  array  radar  as  well  as  commercial  communications 
applications.  Figure  11(b)  shows  measured  two  tone  data  compared  to  measurements  obtained 
with  fREEDA  transient  analysis  and  Agilent  ADS’  harmonic  balance  analysis.  Note  that  in 
transient  analysis  we  can  obtain  higher  dynamic  range  than  obtainable  in  the  industry’s  leading 
commercial  harmonic  balance  simulator.  One  hallmark  of  the  work  that  we  do  is  comparison 
with  measurements  wherever  possible. 
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Output  Power(dB)  at  the  third  order  intermodulation  products 
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Figure  1 1 .  Measured  and  simulated  results  for  an  X-band  MMIC:  (a)  MMIC;  and  (b)  simulation 
results. 
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2.4  Electromagnetic  Development 


2.4.1  Electromagnetic  Aware  Physical  Design  Kit  (EM-PDK) 

EM-PDK  (Java)  is  both  a  stand  alone  program  and  integrated  into  Electric.  Electric  is  a  GNU- 
licensed  IC  design  tool.  We  have  integrated  EM-PDK  into  Electric  and  is  distributed  as  part  of 
Electric.  It  is  written  in  the  Java  language  and  operates  in  two  distinct  modes: 

•  Stand  alone  mode: 

o  Parses  a  layout  file  and  extracts  electromagnetic  (EM)  features  that  are  important 
in  electromagnetic  analysis,  it  is  EM  aware, 
o  Creates  cif  output.  (Metals,  dielectrics  thicknesses  and  properties  only), 
o  Pass  to  electromagnetic  simulator  simulator. 

•  Electric  mode 

o  identifies  local  reference  groups  and  assigns  local 
o  reference  terminals.  Assigns  terminal  numbers, 
o  creates  spice  file  +  em-aware  layout  file 

In  the  Electric  mode  a  tool  flow  using  EM-PDK  is  shown  in  Figure  12. 


Figure  12.  Integrating  fREEDA  and  Electric  usin  EM-PDK. 
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2.4.2  Modifications  to  Electric  Editor 


Electric  Editor  (or  known  simply  as  Electric)  is  a  powerful  layout  tool  that  is  GNU-licensed.  It  is 
widely  used.  It  handles  MOS-Bipolar  transistors,  MOS-Bipolar  Schematics,  and  HDLs.  It  has  a 
common  Layout-Schematic  Interface,  see  Figure  13.  It  supports  various  interface  standards 
including  GIF,  GDS  II,  VHDL,  Verilog.  NCSU  made  a  number  of  modifications  to  include 
naming  of  nodes  at  points  of  connection  between  schematic  and  layout;  support  for  materials 
properties  added  (Layer  Thickness;  Loss  Tangent;  and  Dielectric  Constant)  All  modifications 
have  been  included  in  the  current  Electric  Editor  release. 


r*s»j  tscM  i-iM.  t»» 


Figure  13.  Screen  Capture  of  Electric. 


2.5  IBIS  Macromodeling  Development 

IBIS  is  a  behavioral  macromodeling  language  that  has  been  widely  adopted  by  industry,  see 
http://www.eia.org/eig/ibis/ibis.htm  .  The  structure  of  the  model  is  shown  in  Figure  14.  One 
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reason  why  this  model  has  been  so  successful  is  the  availability  of  an  automatic  macromodel 
generation  tool.  This  tool  is  S2IBIS  (Spice-to-IBIS).  The  development  of  version  3  of  the  tool, 
S2IBIS3,  was  supported  by  NeoCAD.  S2IBIS3  is  being  downloaded  from  the  NCSU  website  at 
the  rate  of  10  per  day  currently.  Brief  highlights  are 

•  Spice  To  IBIS  Macromodel  Converter 

•  Supports  HSpice  &  Spice3 

•  Written  in  Java  (platform  independent)  Tested  Windows,  Unix,  Linux 

•  Provides  heuristic  convergence  assistance. 

•  More  than  1,000  downloads  since  2/1/2004.  (Now  10  per  day) 

•  Version  3  good  for  2  GHz  digital  systems 

•  IBIS  model  in  fREEDA  (first  as  open  source) 

IBIS  is  a  notoriously  difficult  model  to  implement  in  a  simulator  and  there  are  only  a  few 
implementations  in  Spice.  By  providing  an  open  source  implementation  we  have  told  the  world 
how  to  implement  the  model. 

At  the  start  of  NeoCAD  the  start-of-the-art  was  S2IBIS2  the  development  of  which  was  also 
supported  by  DARPA.  S2IBIS2  was  released  1996  and  was  written  in  the  C,  bison,  perl 
languages.  This  resulted  in  platform  portability  problems  and  also  problems  assosicted  with 
operating  system  upgrades.  S2INIS2  could  develop  models  that  were  good  to  300  MHz 

S2IBIS3  is  compliant  with  IBIS3,  version  3  of  the  IBIS  spec.  It  is  written  in  JAVA  and  so  it  is 
platform  independent.  It  was  tested  by  the  EIA/IBIS  group  of  companies  including  mentor 
graphics  and  cadence.  S2IBIS3  supports  Version  3  of  the  IBIS  specification  and  is  good  to  2 
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GHz.  Today  almost  no  printed  circuit  board  system  or  module  can  be  designed  without  IBIS 
models  of  ICs. 


Figure  14.  The  IBIS  model  of  a  driver. 


Nearly  every  semiconductor  vendor  uses  S2IBIS  to  generate  IBIS  models  that  are  generally 
made  available  on  vendor  websites.  There  are  at  least  300  vendors.  Some  of  these  are  shown  in 
Figure  15. 
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Figure  15.  S2IBIS  Semiconductor  vendor  users. 


20 


We  have  received  hundreds  of  testimonials  regarding  S2IBIS  (not  always  in  the  best  English). 
Here  are  a  few: 

•  Hello,  my  name  is  Albert  Kim.  I  am  DRAM  engineer  in  Korea. 

I  would  like  to  say  "Thank  you  for  your  s2ibis3". 

Thanks  to  your  tool,  I  could  save  my  time. 

•  “Great  work  has  been  done,  thank  you  and  good  luck”  -  Christophe  Robichon  (Atmel) 

•  “Toshiba  appreciates  what  you  are  doing!!”  -  Doran  David  (Toshiba) 

•  “Thank  you  for  providing  s2ibis3.”  -  Bob  Weppler  (Rockwell) 


We  also  undertook  studies  to  a  certain  the  accuracy  of  the  IBIS  model  produced.  Figure  16. 
Compares  the  simulated  results  using  Spice  netlist  and  an  IBIS  for  an  off-chip  driver. 


(a) 


(b) 


Figure  16.  Comparison  of  circuit  level  model  and  IBIS  model:  (a)  an  8  inverter  IC  off-chip 
driver;  and  (b)  comparison. 
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2. 6  Advanced  Device  Modeling 


In  fREEDA  many  devices  have  been  modeled  for  the  first  time.  Several  classes  of  models  were 
developed  and  much  work  was  put  into  modeling  types  of  effects  for  the  first  time.  One  first  was 
modeling  the  tunneling  phenomenon  is  an  extremely  nonlinear  process.  A  MOSCAP  is  one  of 
devices  that  was  modeled  and  a  cross-section  view  is  shown  in  Figure  17.  This  captures 
tunneling  through  a  thin  oxide  layer  which  a  particular  problem  in  sub  1 00  nanometer  transistors. 


Figure  17.  Cross-section  of  a  MOSCAP  that  exhibits  tunneling  for  thin  oxide  layers. 

As  is  typical  in  modeling  in  fREEDA  considerable  effort  goes  into  modeling  effects  for  the  first 
time.  Once  a  phenomenon  has  been  captured  it  can  be  replicated.  In  this  case  models  were 
developed  for  Tunnel  Diode,  MOSCAP,  TFT,  and  Molecular  Diode  which  all  exhibit  similar 
tunneling.  This  is  the  first  time  these  devices  have  ever  been  modeled  in  a  simulator. 
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2. 7  Molecular  Device  Modeling 


Organic  molecules  exhibit  a  change  in  their  equilibrium  energy  levels  when  biased  such  that 
current  flows  through  them.  Molecular  device  modeling  is  the  modeling  of  this  transport 
mechanism.  It  involves  measuring  and  simulating  the  conductance  characteristics  of  a  nanopore. 
A  nanopore  consists  of  a  few  hundred  to  a  thousand  molecules  sandwiched  between  two 
electrodes  of  a  few  nanometers  width.  Such  a  device  has  never  been  modeled  in  a  simulator 
before  and  such  modeling  is  essential  in  the  development  of  molecular  electronics.  The 
molecular  device  modeled  is  shown  in  Figure  1 8  and  the  measured  and  modeled  comparisons  are 
shown  in  Figure  19. 


Figure  18.  Molecular  device  model 
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0<perimental  result 

M.  A.  Reed,  C.Zhou,  C.  J,  Muller.  T.  P.  Burgin  and  J.  M.  Tour, 
Conductance  of  a  molecular  junction,"  Science,  Vol.  27B,  pp. 
252-253,  Oct.  1997. 

-  C-V  characteristics 

~ —  i-V  characteristics 


Model  1  (Moiecuie1.cc.  Moiecuiel.h) 
Simulated  result  from  Model  1  in  tREEDA. 


Model  2  (Molecule2.cc.  Molecule2.h) 

Simulated  C-V  and  I-V  result  from  Model  2  in 
/REEDA.  Compared  to  Modei  1,  the 
characteristics  here  capture  the  quantum 
mechanical  nature  of  the  device  more 
accurately. 


Figure  19.  Simulated  results  for  a  molecular  device. 


The  turmel  process  inherent  to  molecular  electronics  cannot  be  expressed  as  a  current-charge- 
voltage  expression.  We  must  use  state  variables  to  implement  tunneling  correctly. 


2.8  Delay  Modeling 


Spice  handles  only  short  (<  4  time  step)  time  delays  using  a  polynomial  interpolation  model. 
fREEDA  has  no  limit  achieved  through  the  implementation  of  an  appropriate  error  measure.  The 
response  an  indefinite  time  in  the  past  is  used  as  a  state  variable.  This  can  only  be  implemented 
using  a  state  variable  based  simulator.  A  TWT  amplifier  was  used  to  validated  fREEDA’s  ability 
to  handle  models  with  long  time  delays  as  shown  in  Figure  20.  This  feature  has  also  been 
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implemented  in  many  transistor  models  in  fREEDA  and  so  we  model  time  delay  in  microwave 
transistors  accurately. 


Time  (ps) 


(/) 


Time  (ps) 

Figure  20.  Demonstration  of  fREEDA’s  ability  to  model  arbitrary  time  delays.  Shown  are  the 
simulated  input  and  out  put  responses  of  an  X-band  MMIC. 


2.9  Details  of  Modeling  of  a  an  X-band  radar  system 


The  final  demonstration  of  the  capabilities  of  the  tools  that  we  developed  were  demonstrated 
using  the  simulation  of  a  radar  system.  A  radar  signal  is  similar  to  that  of  an  FM-chirped  RF 
pulse.  Using  current  tools  it  is  not  possible  to  model  the  nonlinear  dynamics  of  this  system.  The 
design  of  the  RF  section  of  a  radar  currently  utilizes  frequency  domain  tools,  specifically 
harmonic  balance.  A  radar  has  high  power  transient  signals  that  significantly  affect  operation. 
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The  RF  pulse  can  lead  to  spurious  spectral  emissions  and  result  in  transient  heating  effects. 
Capturing  these  effects  in  a  transient  simulation  is  not  possible  with  current  tools. 

The  attributes  required  to  model  a  radar  system  with  compromise  are  shown  in  Figure  21.  At  the 
start  of  the  program  there  was  limited  ability  to  model  the  dynamic  range  of  the  radar  system.  An 
exception  was  using  harmonic  balance  as  in  the  ADS  simulator  but  this  simulator  is  not  capable 
of  capturing  true  radar  effects.  Since  the  program  started  both  Spectre  and  ADS  (with  its  inbuilt 
transient  simulator)  have  improved  their  dynamic  range  in  transient  simulation.  However  they 
do  not  have  the  dynamic  range  of  Spice.  It  is  not  clear  exactly  what  their  essential  dynamic 
range  as  it  appears  that  it  is  limited  by  their  internal  Fourier  transform  algorithm.  For  the 
fREEDA  results  presented  here  we  used  a  special  Fourier  analysis.  With  fREEDA  the  dynamic 
range  is  160  dB  (to  an  accuracy  of  0.5  dB).  If  we  applied  the  same  criteria  to  ADS  and  Spectre 
we  could  say  that  the  range  is  limited  to  120  dB.  The  attributes  required  all  of  the  fREEDA 
developments  mentioned  previously. 

As  fREEDA  is  public  domain  any  company  can  examine  the  code  and  determine  how  to  achieve  the 
same. 

Feature  Initial  State  of  The  Goal  Achieved 

Art 

Dynamic  Range  40-60  dB  SPICE  140  dB  160  dB 

>  120  dB  ADS 

Multi  Physics  Limited  Thermal  /  EM  /  Circuit  Thermal  /  EM  /  Circuit 
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Time  Delay 


2  or  3  time  steps  Unlimited 


Unlimited 


Figure  21.  Initial  State  of  the  Art,  Goal  and  Achievement  of  fREEDA  as  used  I  transient  radar 
system  simulation. 

While  were  unable  to  obtain  models  for  all  of  the  components  of  a  radar  system  so  we 
constructed  one  using  available  models  of  actual  devices  and  idealized  models.  The  radar 
system  model  considered  is  shown  in  Figure  22.  It  consists  of  a  Linear  FM  Chirp  source  pulse 
modulated  with  center  frequency  10  GHz  and  lOOMHz  bandwidth.  It  included  the  X-band 
Filtronic  Solid  State  MMIC  LNA,  LMA  41 1  with  18dB  gain  (considered  previously).  The 
MMIC  was  followed  by  a  Traveling  Wave  Tube  (TWT).  A  circuit  level  model  of  the  MMIC 
LNA  was  used  and  as  seen  previously  simulations  accurately  capture  measured  data.  The 
Bandpass  Filter  removes  harmonic  content  from  the  MMIC  output.  Transient  analysis  of 
bandpass  filters,  especially  narrow  band  bandpass  filters  is  problematic.  We  have  developed  a 
technology  based  on  Z-domain  techniques  that  enables  any  order  filter  to  be  modeled  in  transient 
analysis.  The  TWT  amplifier  was  model  by  a  behavioral  model  using  a  novel  behavioral  model 
and  based  on  measurements.  One  characteristic  of  the  TWT  amplifier  is  the  very  long  time 
delays  involved.  Modeling  this  structure  in  fREEDA  amounts  to  the  first  time  a  TWT  has  ever 
been  modeled  in  a  transient  simulator.  There  are  many  firsts  here  and  we  will  write  journal 
papers  up  to  cover  all  of  these  advances.  Our  primary  objective  recently  has  been  to  achieve  the 
final  phase  milestones. 
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Various  results  are  presented  in  Figures  23  through  29  for  the  transient  waveforms  and  the 
spectra.  The  spectra  observed  could  only  be  captured  by  first  doing  transient  analysis.  In  Figure 
22(b)  missing  features  are  indicated. 
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Radar  Xmit  System  Block  Diagram 


Chirp  Source:  FM  Linear  downchirp  from  1.5  to  0.5  GHz 
LO  Source:  10  GHz  sine  wave 
Summer  and  Square-Law  mixer  -  both  behavioral 
BPF:  Z-domain  cascade  behavioral.  -1  dB  flat  passband 
from  8.5  to  9.5  GHz,  -  40  dS  stopband  edges  at  8.35  and 
9.75  GHz. 

Driver:  Filtronlc  LMA41 1  MMIC  -18  dB  gain  (circuit  model) 
PA:  HP495A  TWTA  -  20  dB  gain  (behavioral) 

Output  Signal:  Linear  FM  Upchirped  from  8.5-9.5  GHz. 


Still  To  Do... 

White  Noise  source  (based  on  gaussian  random 
variable) 

■  User-specified 

•  To  facilitate  receiver  modeling 

•  Presently  generating  random  variates,  but  not  at 
correct  power  spectral  density 

Channel  Model 

•  Time  delay  plus  insertion  point  for  white  noise  source 
model,  to  facilitate  receiver  modeling 

New  -  Chebychev  Z-domain  filter  model 

•  Based  upon  Butterworth  model 

■  Turnaround  time  to  develop  should  be  short 


(a) 


(b) 


Figure  22.  Canonical  radar  system:  (a)  model;  and  (b)  features  missing. 


Linear  FM  Downchirp  Source  (1)  Linear  FM  Downchirp  Source  (1) 

From  0  to  nt.  PiJierepetitioa  period  ii  15  n>,  chiry  i(  10  nc 

HIT  fU  lI*P  out 


(a) 


(b) 


Figure  23.  Waveforms  at  node  (1):  (a)  short  time;  and  (b)  long  time. 
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Sinusoidal  Local  Oscillator  (2) 

Frcrm  0  to  100  pt,  i.e  two  periods  oflocjd  otriliiilor. 

Din  out 


Sinusoidal  Local  Oscillator  (2) 


ritfv_2  irm  cut 


0  W<TO 


(b) 


Figure  24.  Waveforms  at  node(2):  (a)  short  time;  and  (b)  long  time. 


Summed  Chirp  and  LO  (3)  Summed  Chirp  and  LO  (3) 


IVom  0  to  16  ns.  Kate  thxt  chirp  pertirbKles  the  cxm'cr  envelope 
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(a)  (b) 


FFT  of  Summed  Chirp  and  LO  (3) 
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Figure  25.  Waveforms  at  node(3);  (a)  short  time;  (b)  long  time;  and  (c)  spectrum. 
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FFT  of  Square-Law  Output  (4) 
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Figure  26.  Waveforms  at  node(4):  (a)  short  time;  (b)  long  time;  and  (c)  spectrum. 
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Figure  27.  Waveforms  at  node  (5):  (a)  short  time;  (b)  long  time;  and  (c)  spectrum. 
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Figure  28,  Waveforms  at  node  (6):  (a)  short  time;  (b)  long  time;  and  (c)  spectrum. 
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Figure  29.  Waveforms  at  node  (7):  (a)  short  time;  (b)  long  time;  and  (c)  spectrum. 

We  are  at  the  point  where  we  can  model  any  radar  system  at  any  level  of  accuracy.  Our  one  big 
remaining  problem  is  the  implementation  of  a  BSIM  MOSFET  model  which  is  essential  if  we 
are  to  model  the  types  of  radar  systems  that  Wright  Patterson  are  interested  in.  Implementing  a 
BSIM  model  was  not  in  the  initial  scope  of  this  project  and  has  proven  to  be  problematic.  We 
are  secured  corporate  funds  to  continue  this  part  of  the  work. 
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3  Multigrid-Enhanced  EM  Modeling  of  the  Power  Grid  On  Chip 


This  EM  modeling  based  on  Maxwell  equation  using  a  finite  volume  model.  The  key  attributes 
are: 

•  Model  developed  directly  from  physical  structure 

•  Cumbersome  and  error-prone  extraction  of  [L]  and  [C]  avoided 

•  Rigorous  modeling  of  electromagnetic  effects 

•  In  addition  to  power  switching  noise  analysis  it  enables  prediction  of  power  grid-induced 
EMI  between  different  blocks  on  the  chip 

The  On-Chip  Power  Grid  Transient  Simulator  combines: 

•  Electromagnetic  rigor 

•  Comprehensive  modeling 

•  Includes  impact  of  semiconductor  substrate 

•  SPICE  Compatible.  Initially  we  thought  that  we  could  only  do  this  in  fREEDA  and 
although  it  is  most  efficiently  done  there  we  developed  a  scheme  so  tat  it  can  be 
implemented  in  Spice. 

•  Direct  implementation  of  SPICE  models  for  non-linear  drivers,  decoupling  caps,  etc. . . 

•  Supports  convenient  interfacing  with  models  for  the  off-chip  power  distribution  network 

•  Supports  both  simultaneous  switching  noise  and  power  grid-induced  interference  prediction 

Results  are  documented  in  Figures  30  through  39. 
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FV-Based  Modeling  Methodology 

Finite-volume  discretization  of  Maxwell's  equations 


•  Micron-size  cross-sectional  dimensions  and  regular 
layout  of  the  grid  exploited  to  contain  model  complexity 

-  Grid  size  of  the  order  of  grid  feature  size 

-  Assignment  of  unknown  electric  &  magnetic  fields  in  space 
dictated  by  the  electromagnetic  effects  that  must  be  captured 
for  accurate  simulation 

•  Ohmic  loss  in  the  wires 

•  Inductive  effects  during  switching 

•  Capacitive  coupling  and  common  impedance  coupling  For 
grid-induced  interference  calculation 


Figure  30  FM-based  modeling  methodology. 

Implementation  of  Variable-Size  Grid 


The  choice  of  grid  coarseness  is 
dependent  on  the  simulation  objective  and 
the  desired  accuracy. 

-  Switching  noise  simulation  only;  Coarse  grid 

-  Power  grid-induced  interference:  Finer  grid 


Figure  31.  Implementation  of  variable-size  grid. 


Compatibility  with  SPICE 

The  state-space  form  of  the  discrete  model, 
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is  of  the  same  form  with  the  MNA  formalism  used  in  SPICE 

Hence,  it  facilitates  direct  incorporation  of  lumped  circuits 
and  behavioral  models  for  drivers  and  receivers 

•  e ,  h  :  the  discrete  unknown  fields 

•  G.  R,  L,  C.  Dp,.  Dg :  sparse  matrices  (dependent  on  material  and 
geometric  properties  of  the  structure) 

•  \q,  Vg ;  voltage  &  current  sources  connected  to  the  grid 


Figure  32.  Compatibility  with  Spice. 

“Cartoon”  of  the  On-Chip 
Grid 
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Figure  33.  On-chip  grid. 
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Figure  34.  Metal- 1  layer. 

Additional  Features  of  Solver 

•  Convenient  interface  to  models  for  the  off- 

chip  power  grid 

-  Can  be  incorporated  either  in  terms  of  SPICE 
net  lists  or  in  terms  of  a  matrix  rational 
function  representation  of  their  multi-port  form 

•  Modeling  of  the  semiconductor  substrate 

-  Effected  through  position-  and  frequency- 
dependent  surface  impedance  relationships 
cast  in  terms  of  rational  functions 


Figure  F.  Features  of  solver. 
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Impact  of  semiconductor  substrate 
on  power  grid  switching  response 


Case  0:  Semiconductor  substrate  modeled  as  PEC 

Case  1 :  t,  =  198  /jm,  a,  =  104  S/m,  t2  =  2  fim,  and  Oj  =  10  S/m 

Case  2:  t,  =  200  jim,  a,  =104  S/m,  tj  =  0  {no  epi) 


Figure  36.  Substrate  impact. 


Impact  of  semiconductor  substrate.. 


Figure  37.  Substrate  impact. 


Transient  EM  Modeling  of 
Power  Switching 

Visualization  of  on-chip  supply  voltage  disturbance 
during  simultaneous  switching  at  all  nodes 

t;.  simiAanaous  switching 


Figure  38.  Model  of  power  switching 
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Figure  39.  Model  of  substrate  noise. 


4  Optb-ElectronicModeling  using  fREEDA 


4.1  Introduction 


Progress  in  the  field  of  opto-electronic  device  and  optical  component  technol¬ 
ogy  lias  led  everyone  to  believe  that  optics  can,  practically,  meet  the  bandwidth 
demand  in  future.  Efforts  are  on  to  make  all  optical  systems  where  most  of  the 
tasks  e.g.  switching,  routing  and  amplification  are  done  in  optical  domain.  In 
spite  of  all  these  developments  w^e  still  need  drivers  to  turn  on  a  laser  and  so 
we  have  a  lot  of  electro-optical-thermal  interactions  taking  place  right  at  the 
point  where  light  is  actually  generated.  Further,  when  optical  components  and 
devices  are  put  together  in  a  system  prolilems  can  arise  because  of  the  disper¬ 
sive  elements  and  feedback  effects.  So,  before  start  fabricating  the  system  one 
must  study  the  electro- tlicrmal-optical  interactions  and  other  optical  effects 
such  as  feedback  and  effect  of  dispersion.  To  study  such  complex  interactions 
one  need  a  tool  that  can  be  used  to  simulate  the  system  directly  in  terms  of 
equations  representing  these  interactions.  Existing  tools  for  OE  system  simu¬ 
lation  involve  either  solving  Maxwell  eciuations  using  for  example,  the  FDTD 
method  that  is  computationally  intensive  or  use  equivalent  circuits.  Imple¬ 
mentation  of  device  models  using  an  equivalent  circuit  approach  has  several 
drawbacks. 
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We  have  developed  a  new  tool,  fR,EEDA,  capable  of  performing  mnlti-physics 
simulations.  fREEDA  facilitates  OE  device  modeling  with  sets  of  differential 
equations.  This  report  is  organized  as  follows.  In  first  section,  we  discuss  the 
validation  of  fREEDA  for  optoelectronic  simulation  using  experiments  and  via 
comparison  with  other  existing  tools.  Section  II  discusses  the  multimode  VC- 
SEL  modeling  using  fREEDA.  Next,  we  present  the  simulation  results  for  free 
space  optical  system  using  fREEDA.  Optical  feedljack  effects  and  there  im¬ 
pact  on  the  system  performance  ai'e  discusses  in  the  following  section.  Finally, 
we  conclude  with  a  summary  of  results  and  the  future  work. 

4.2  fREEDA  Validation 


In  this  section,  we  discuss  the  vali<lation  of  fREEDA  for  Vertical  Cavity  Surface 
Emitting  Laser  (VCSEL)  modeling  and  simulation  using  the  experiments  and 
CFDRC  physical  VCSEL  model.  We  also  present  the  comparison  of  IREEDA 
results  with  the  existing  tools  such  as  SPICE  for  package  simulations. 

We  chose  VCSEL  for  fREEDA  validation  as  its  optical  properties  depend 
strongly  on  the  ojierating  electrical  and  thermal  conditi(.)ns  [2].  So,  it  makes  tlie 
most  suitable  OE  device  for  validating  the  multi-physics  nature  of  IKEEDA. 
Validation  of  IREEDA  with  the  experiments  an<l  other  simulators  have  been 
done  in  past[9,10].  Wo  u.se  the  VCSEL  model  of  Mena  in  our  simulations  and 
use  the  following  rate  eejuations  [2]: 

^  _  <^(^^0  -  N,)S  _  h 

dt  q  (1  -f  eS)  q 

dS  S  ,dNo  .  G{No-Nt)S 
dt  Tp  T„  (1  -f  eS) 

where  No  is  the  carrier  number,  S  is  the  photon  number,  yV,  is  tra.nsi)arency 
carrier  density,  77,-  is  the  current  injection  efficiency,  /5  is  the  spontaneous  emis- 
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sioii  factor,  e  is  the  gain  saturation  factor  and  //  is  the  leakage  current.  Carrier 
lifetime  and  phot(3n  lifetime  are  t„,  and  Tp  respectively. 

We  model  the  instantaneous  wavelength  of  the  VCSEL  using  the  equation 
given  in  [1]. 


A(t)  -  A„(l 


p{No-Nt)  dn{r-To) 
n  dT  n  ^ 


(3) 


The  instantaneous  wavelength  depends  on  carrier  number  and  temperature. 
To  is  the  ambient  temperature,  p  is  the  refractive  index  change  per  unit  carrier 
number  change  and  n  is  the  refractive  index  of  the  medium.  Fig.  1  shows  a 
comparLson  of  the  experimental  and  simulated  L-I  charticteristics  for  a  sin¬ 
gle  mode  Honeywell  VCSEL  SV3639.  Simulated  results  agree  well  with  the 
experimental  results  over  a  range  of  temperature. 


Fig.  1.  Experimental  and  simulated  L-I  curves  for  VCSEL  at  different  temperatures. 

We  also  used  the  L-I  and  I-V  cliaracteristics  from  the  CFDRC  physical  VCSEL 
model  to  extract  the  VCSEL  parameteis.  Fig  2.  shows  the  CFDRC  VCSEL 
model  along  with  the  parameter  extraction  procedure.  Simulated  L-I  and  I-V 
characteristics  are  compared  with  the  CFDRC  VCSEL  L-I  and  I-V  cliaracter- 
istics  in  fig.  3a  and  fig.  3b  respectively.  There  is  a  good  agreement  lietween 
the  simulated  and  the  CFDRC  pliysical  VCSEL  characteristics. 
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CFnRCVCiiri  CJ1>RC VCSIO. 


Fi^^.  2.  C’l'DRC  VC’SEL  [jliysical  model  niul  parniiieter  extraction  procedure. 


(a)  (b) 

Fig.  3.  Coinpnri.son  of  the  C'l*’DTtC'  VCSEL  cliaracteristics  with  the  ITIEEDA  results 
using  extracted  paratueRns  .  (a)  L-l  characteristics  at  difrereut  temperatures  .  (Ir) 
l-V  characteristics  at  {liHeri'nt  ten iperat art's. 


Wv  also  did  a  .study  of  flio  eh’etrical  interaetioii.s  that  take  [)lac(‘  in  tht'  package.s 
eomiectiug  the  VCSEL  array  to  tli('  VCSEL  driver  array  [Ij.  fREEDA  i>aekage 
results  are  eomi>arefi  with  the  SPICE  results  in  figures  4..'3  and  6.  Fig.  4  aiifl 
5  show  a  eomi>arisou  of  output  optical  pcjwer  fr<»ui  fREEDA  and  SPICE  for  a 
tliin-film  i>aekage  at  IGBps  and  2.5  GBps  respectively.  A  eoinparison  of  the 
output  optical  [xjwer  for  a  flip-chip  i>ackage  at  5GBps  is  presented  in  fig.  C. 
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X  lO"* 


Fig.  4.  Conipari.son  of  fREEDA  and  SPICE  results  for  thin-filni  paclcage  at  iGBps. 


X  iO"’ 


Fig.  5.  Comparison  of  fREEDA  and  SPICE  results  for  thin-film  package  at  2.5GBps. 
Details  of  the  packages  used  in  simulations  are  given  in  [1]. 


4.3  Multimode  VCSEL  Modeling 


Next  we  pre.sent  the  siinulation  results  for  a  multimode  VCSEL.  We  modeled 
the  multimode  VCSEL  using  the  model  described  in  [2].  Following  ecpiations 
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Fig.  6.  Comparison  of  fREEDA  and  SPICE  rcsnlts  for  flip-chip  package  at  5GBps. 


were  used  in  our  simulations: 


dNp  _  Tjjl  ^  ^  ^  Gjjko^o  —  Ihi^i  —  'ykoNt)Sk  _  h  / 
dt  q  r„  ^  (1  +  ErnGnArS’,,,.)  q 

d'Nj  _  -1-  ^  _L  G((pjicoNo  —  (pjkjNj  d^jkoNf^Sk 

Q  Z^TYi  ^'nxk^rt^} 


dt 


Sk  .  Sk 


Tpk 


+  ^ — -[^oA^o  —  ^  biNi]  -I-  y~^ 


T, 


G(XicoNo  —  Ei  ^kiN-i  —  XkoNt)Sk 

(1  +  E  ■ni  ^7nk‘^  m) 


(6^ 


These  siiatially  independent  rate  ecpiations  arc  derived  and  explained  in  [2]. 
Though,  these  equations  are  spatially  independent  they  very  v'cll  capture  the 
spatial  effects,  sucli  as  spatial  hole  l)urning,  required  for  multimode  ox^eration 
through  A^i-  We  rei^resent  the  DC  and  transient  results  for  a  two  mode  VCSEL. 
Fig.  7  shows  the  L-I  characteristics  for  a  two  mode  VCSEL.  From  fig.  7  it  is 
clear  that  the  threshold  for  mode  2  is  higher  than  that  for  mode  1  since  the 
emission  cjf  mode  2  require  a  spatial  hole  to  burn  which  occurs  at  higher 
current.  Once  a.  sx^atial  hole  is  bui’iied,  the  carrier  x^iofile  changes  to  siij^xiort 
the  emission  of  LPn  mode.  A  kink  appear  in  the  L-I  characteristics  of  LPqi 
mode  at  the  i^oint  where  LPn  mode  begins  to  lase,  changing  the  slox^e  of  the 
L-I  curve  for  LPpi  mode. 
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Fig.  7.  L-I  characteri.stics  for  a  multimode  VCSEL.  Continuous  line  shows  LPoi 
mode  and  dashed  line  shows  LPn  mode. 

Transient  results  for  a  two  mode  VCSEL  are  shown  in  figui-e  8.  Here  also  we 
see  that  the  mode  2  start  lasing  at  a  later  instant,  showing  the  impoi'tance 
of  SHB  in  the  emission  of  LPn  mode.  Once  tlit;  output  of  mode  1  goes  down 
it  comes  up  again  which  is  .seen  as  a,  little  hump  in  the  LPoi  power  output. 
These  secondary  pulsations  occur  due  to  the  diffusion  of  carriers  towards  the 
spfitial  hole  once  the  device  is  switched  off.  These  results  show  the  capability 
of  fREEDA  in  simulating  the  complex  spatial  behavior  of  a  multimode  VCSEL 
using  the  system  of  ecpiations.  To  simulate  the  similar  behavior  in  other  circuit 
simulators  one  need  to  generate  an  equivalent  circuit  corresponding  to  these 
equations.  Further,  different  circuit  models  for  different  types  of  analysis,  such 
as  DC,  transient  and  AC,  are  required  in  other  circuit  simulators,  wherea.s  any 
type  of  ana.h^sis  can  be  performed  in  fREEDA  using  the  same  set  of  equations. 


4.4  Free  Space  Propagation 


This  section  discusses  the  modeling  of  a  free  space  optical  system  using  fREEDA. 
We  model  the  free  space  propagation  using  the  gaussian  beam  propagation 
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4.5, 


Fig.  8.  Transient  characteristics  for  a  multimode  VCSEL.  Continuous  line  shows 
LPqi  mode  and  dashed  line  shows  LP\i  mode. 

model.  For  gaussian  proi^agation  the  beam  width  at  a  distance  z  is  given  by: 

w{z)  =  woy/l  +  {zjzrY  (7) 


where  wo  is  tlie  beam  waist  size,  Zr  is  the  rayleigh  range.  Rayleigh  range  Zt  is 
related  to  the  wavelength  A  and  the  beam  waist  Wo  by  the  following  ecpiation: 


r  =  TTte^/A 


(8) 


Fig.  9  shows  a  simple  free  sjrace  optical  system  imaging  an  inpr.it  beam  to 
output  Vreain  using  a  lens. 


Fig.  9.  System  diagram  for  free  space  gaussian  propagation. 


However,  finite  size  of  the  lenses  and  other  apertures  clip  the  beam  leading 
to  power  reduction  and  diffraction  effects[3].  Power  reduction  caused  by  the 
finite  circular  aperture  can  be  modeled  using  ecpiation: 


-Priem  =  P(1  -  exp(-2P)) 


(9) 
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where  Ic  =  is  t  he  rat  itJ  of  the  apert  ure  diameter  to  tlie  spot  size  at; 

tlie  aperture  input.  Presenee  of  the  aperture  modifies  the  beam  widtii  to: 

=  »’o(l  -  exp(-2A-^).eo*'{/)A-'-^/(27r)))  (10) 

aiui  tlie  new  ra\-leigli  range  to: 

Z,n-  =  7TW:,„j„/X  (11) 

where  p  =  27r(.:/ii,-),  z,-  is  the  initial  rayleigh  range  and  z  rejiresent  tlie  distance 
from  waist  to  tlie  aiiertnre.  We  imaged  a  lasfw  lieani  of  sjiot  .size  iCo  =  0/nn 
using  a  leii.s  of  focal  length  f  =  50//m.  The  laser  was  plated  at  lOO/zm  and  the 
spot  size  and  i^ower  was  detectefl  at  100/mi.  To  observe  the  effect  of  clipping 
we  varied  the  lens  diameter  from  I0//m  to  33/nn.  Fig.  10  and  11.  show  the 
effect  of  cliiiiiing  on  fletected  spot  size  and  power.  Detected  power  aiul  si>ot 
size  did  not  show  any  change  for  an  aperture  diameter  of  33/nn.  But  a.s  the 
aperture  size  is  reduced  further,  spot  size  increa.ses  due  to  diffraction  leatling 
to  reduced  deti'cted  iiowcr. 

D  D  □  □ 

dims=3L3um  dlHti*=22ttni  dlBis^l5uin  illitisr  lOiim 

Fig.  10.  FfTect  of  lots  ai>erlure  size  on  spot  diagram  at  cleteetor  plane. 

We  also  modeled  the  free  space  propagation  using  the  following  ray  tracing 
etjuations; 


Q  —  f{coft(U)  —  usin{lJ) 

{12} 

(13) 

sh}{fi)  =  ns//)  (/)/«! 

m 

U,  =  i;  _  /  +  /, 

(15) 

49 


Fig.  11.  jjowt'r  for  Ions  npcrtiirr  sizos  showing  rt'diU'Tioii  in  dr^ 

toftod  power  ns  tlie  lens  npertiiie  is  roduef'd, 

(h  =  Q{ros{lh)  +  ros{h))/{ros{V)  +  r-o.s'(/))  (Ifi) 

\vher<«  Q  is  tlw'  (lisfnnee  from  the  vertex  of  the  surface  to  the  hu  icleiit  jay  and 
perpi'iidieular  to  the  ray.  Qi  is  tlie  tlistaiiee  from  tlu'  vert('x  of  the'  surface  to 
the  refraeted  ray  and  pci  jieudieuiar  to  tlie  ray,  H  is  the  input  ray  height,  U 
is  the  ininit  ray  angle,  I  is  the  angle  of  iiu-idenee.  /i  angle  of  refraetion  and  u 
objt'et  distauee. 

Ray  tracing  is  good  foi-  the  fiist  onler  optical  design  a.s  it  IningK  out  the 
aberrations  present  in  tin'  systcmi.  Figures  12n  and  12h  show  tlie  c^tfeet  of 
plano-convex  Icms  on  the  spheric  al  ahc'iration.  We  sc'e  that  sjiherieal  ahen  ration 
is  more  win'll  plane  snrfac-e  is  placc'cl  towards  tlu'  .source,  which  is  tistially  tlu' 
ease'.  F^ie.si^iu'e  of  field  c  urvature  is  .sliown  in  lig.  13. 


Thc’  effects  (jf  aherratious  and  diffraetioii  c'an  lie  reduced  at  an  early  stage  of 
ch'.sign  by  choosing  the  lenses  corrected  for  aberration  anrl  witli  large  njiertnrc'. 
But  prohlc'ins  can  still  arise  diu'  to  the  prc'.sc’ncc'  of  dispersive  ek'iuc'nts.  W'hc'ti 
a  laser  is  dircu  tly  modulatcnl  it  givc's  rise'  to  the  ringing  noise*  and  wavek'ngtli 
cliirp  [1].  Fig.  14  shows  a  system  having  a  diftraction  grating  as  a  disjier- 
sive  ek'mc'iit .  Bec  ause  of  the  wavc’lc'iigtli  chirp,  the  grating  eaiisc'S  a  temporal 
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(h)  {!>) 

Fig.  12.  Rjiy  traco  tliagram  lor  plano-convex  lens,  (a)  With  convex  -smfacc  toward.s 
the  source,  (h)  With  piano  surface  tow'aixls  the  source. 


'll  «•  i  i!a  »  14  u  I*  a 

*  4 

a  >■* 


Fig.  l.'h  Ray  trace  for  a  l>inconvex  len-s  showing  field  cnrvanire. 

\ariati(»u  in  file  beam  ])«:>sition  on  the  deleetor  plane  as  shown  in  fig  15. 

Tliis  temporal  variation  in  the  position  on  the  tletector  plane  git^cs  rise  to 
d(?te(  lor  crosstalk  in  a  fret!  space  optical  itit(Tconnect(FS01).  Fig.  IG  shows 
the  sptjf  size  on  tliffereut  detector  po.sitioii.s  as  a  function  of  time.  Optical 


Datfietor 

Array 


Fig,  14.  Sysk'in  (lingraiii  for  givillTig  slK'ctroiiit'tcr. 


Fig.  15.  Variat  ion  of  spot  |>ositioii  on  rlnt.ootor  «irra,v  plant'. 
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Fig.  16.  Spot  diagram  for  variation  of  si>ot  po,sition  on  tletcctor  array  witli  time, 
power  iiioasmod  Ijv  tlio.se  dotootor.s  is  shown  in  lig.  17. 

Howt'Vtu’.  wavelength  chirp  can  ho  oliminatod  using  external  niodulatioii.  Fig. 
18  shows  an  electro-optic  modulator,  whi'ro  the  aj>[)licntion  of  externa lly  ap- 
plit'd  voltagt'  caust'.s  a  pha.se  shift  between  the  field  in  two  arms  of  a  Maeh- 


52 


350| 


Fig.  17.  Output  optical  j>o\vcr  v.s.  time  for  «liH'eretit  fletectors  in  the  cleteetor  Hrray. 


Fig.  IH.  A  .simple  EU-etrtj-optie  niodulntor  fur  external  moilnlntion. 

Zeliiider  iiitcrferouietor.  Tlic  output  of  a  Macli-Zeluidcr  modulator  cau  be 
oxjjres.sed  in  terms  of  the  applied  voltage  using  the  following  crpiation: 

/,>,e  =  /,„(r<^s(7rV7(2  V;)))2  (17) 

where  I„„,  is  the  (jutput  intensity.  /,„  is  the  input  intensity,  14  the  wltnge 
required  for  tt  phast'  shift  and  V  is  tlie  applied  voltage. 

Fig.  10a  and  19b  shows  the  output  optical  power  and  wavelength  with  and 
without  the  external  uiofhdalion.  Optical  power  and  the  wavc'length  for  the 
direct  nuxlnlation  case  .sliows  ringing  and  chirp  but  the  extc'rnal  modulation 
reduces  it. 
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—  —  External  moduiation  I 
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(a)  (b) 

Fig.  19.  Comparison  of  external  and  direct  modulation,  (a)  Output  optical  for  direct 
and  external  modulation,  (b)  W^vvelengtli  vs.  time  for  direct  and  external  rnoduUv 
tion. 

The  impact  of  external  modulation  is  more  clear  from  tlie  eye  diagrams  shown 
in  fig.  20a  and  201).  Fig.  20a  shows  that  the  eye  opening  is  less  in  case  of  direct 
mochilation  due  to  the  ringing  noise,  whereas,  the  external  modulation  shows  a 
good  eye  opening  in  fig.  201).  Further,  use  of  the  external  modulation  eliminate 
the  detector  crosstalk  that  occurs  in  the  pre,sence  of  chirp.  Spot  position  on 
the  detector  plane  is  compared  for  the  direct  and  external  modulation  in  fig. 
21  as  a  function  of  time.  For  external  modulation,  the  spot  position  does  not 
change  with  time,  whereas,  for  the  direct  modulation  spot  position  jumi)s  from 
one  detector  to  another  with  time. 

Thus,  fREEDA  enables  modeling  and  design  of  oi)tical  systems  using  the  sim¬ 
ple  set  of  equations.  To  perform  a  similar  system  study  using  other  circuit 
simulators,  recpiire  converting  each  and  every  ecpiation  into  an  equivalent  cir¬ 
cuit  [4] .  We  see  that  one  can  study  the  various  interactions  that  take  place  in 
a  system  using  fREEDA  but  the  most  unicpie  feature  of  fREEDA  is  its  ability 
to  make  the  optical  feedback  study  possible.  fREEDA  has  a  delay  element 
that  makes  it  possible  to  take  into  account  the  reflection  from  each  optical 
surface  bac'k  to  the  laser  source.  Next  section  is  devoted  to  optical  feedback 


54 


(»>)  {10 

Fi}'.  20.  CV>nipari.son  of  t'Vf'diusroins  for  external  and  dinx-t  modidalioti.  (a)  FjVc 
diagram  for  dirwd  modulatioti  shows  les.s  eye  opening,  (h)  E^'e  diagram  for  external 
modulation  show  good  cyc'  opening. 
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Fig.  21.  CVjmparisou  of  .spot  diagrams  on  detector  plane  for  direct  and  external 

modulation. 

study  using  fREEDA. 

4.5  Optical  Feedback  Study 

Optical  feedback  plays  an  important  role  in  optic:al  communication  systems. 
Even  a  small  amount  of  feedback  can  cause  instability  in  the  laser  output 
resulting  in  system  performance  degradation.  But,  one  can  also  use  optical 
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feedback  for  useful  applications  e.g.  strong  feedback  can  give  rise  to  single 
mode  operation  and  line  narrowing.  So,  one  should  consider  reflection  from 
optical  surfaces  while  designing  an  optical  system.  We  model  the  feedback 
using  Lang-Kobayashi  ecpiations  [5]. 
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Where  r  =  is  tiie  round  trip  time  of  external  cavity,  k  —  \J Bext  / R, 

Rext  is  tlie  reflectivity  of  external  reflector,  R  is  tlie  reflectivity  of  laser  mirror 
and  T£,  is  the  round  trip  time  of  internal  cavity.  Optical  feedback  has  been 
studied  extensivedy  by  various  groups  and  five  regimes  of  feedback  have  been 
identified.  Fig.  22  show  various  regimes  of  optical  fecidl)ack.  Fig.  22  lias  the  out¬ 
put  optical  power  and  the  corresponding  phase  plots.  We  see  that  the  output 
optical  power  aiiiiroaehes  stable  point  for  low  reflectivity  and  sliort  distance. 
As  the  distance  and  reflectivity  is  increased,  output  power  exhibit  a  limit  cycle 
followed  by  chaotic  output. 

To  our  knowledge  none  of  tlu'  groups  have  studied  the  impact  of  the  ojitical 
feedback  on  system  performance.  Here,  we  discuss  the  impact  of  the  optical 
feedback  on  system  performance.  Fig.  23  shows  an  imaging  system,  with  ojiti- 
cal  feedback  from  different  reflective  surfaces.  We  calculate  the  resultant  field 
at  each  surface  and  then  projaagate  it  back  to  the  laser. 


In  our  study,  we  c;onsider  the  reflection  from  single  lens  and  tw'o  lens  system. 
The  outiiut  optical  power  and  output  wavelength  for  these  two  systems  are 
sliown  in  fig.  24  and  fig.  25. 
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Fig.  22.  DiftV'ifiit  rcj^imos  of  opticnl  ft*(xll)ack. 


Fin. '23.  Two  Icni.s  sy.stcnu  tlopictiun  optkfil  fo<Hlh}u  k  from  tUfferoiit  comjtonents. 

Wo  soc  that  the  optical  i)owor  iiiul  waA'cleiiRth  .show  iii.stal>ilitj'  cotilirming 
that  even  a  small  aniomit  of  ftx‘(l!)ack  can  (lenrade  the  system  perfornumce. 
In  previcHjs  section.  \vc  saw  that  by  using  the  external  modniation  one  can 
eliminate  tlie  cliirp  and  imjjiove  eye  opening  l)nt  there  wo  didirt  consider  the 
eff(‘ct  of  fecN  iback.  To  sec*  how  optic  al  fcH'ilback  imijact  the  system  performance, 
we  ecinsider  the  grating  spectrometer  system  shown  in  fig.  14  and  lookt'd  at 
the  outj>nt  optical  power  on  a  fixed  detector  in  tlie  detector  array  as  sliown  in 
fig.  26. 
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Fig.  24.  Effect,  of  ojotical  feedback  cjn  output  optical  power  in  a  .single  lens  and  two 
lens  system. 
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Fig.  25.  Effect  of  optical  foodback  on  output  wavelength  in  a  single  lens  and  two 
lens  system. 

The  output  optical  powder  varies  randomly  with  time  on  the  detector  because  of 
instability  in  wavelength.  If  there  had  been  no  feedback  and  we  were  using  the 
external  modulation,  then  tlie  output  spot  should  have  stayed  at  one  detector 
as  shown  in  fig  21.  Pig.  26  shows  the  output  optical  power  on  a  fixed  detector 
for  no  feedback,  feedback  from  a  single  lens  and  feedback  from  two  lenses. 


Further,  we  noticed  that  for  small  feedback  and  small  distance,  within  coher¬ 
ence  length  of  laser  output,  outinit  ]>ower  shows  a  sinusoidal  behavior  with 
period  ecjual  to  half  the  wavelength.  We  also  noticed  that  the  feedback  causes 
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Fig.  26.  Effect  of  feedback  on  detector  output  for  single  lens  and  two  lens  system. 

a  change  in  the  teriniiial  voltage.  Fig.  27  shows  a  plot  of  voltage  vs.  reflectivity. 
The  terminal  voltage  reduces  as  we  increase  the  reflectivity  because  increased 
reflectivity  causes  more  charge  carriers  resulting  in  reduced  input  impedance 
and  thus,  the  voltage.  One  can  utilize  this  variation  in  voltage  with  feedback 
for  making  comijact  imaging  or  readout  systems  [7,8].  The  cliarige  in  the  in¬ 
put  impedance  with  optical  feedback  can  cause  impedance  mismatch  between 
the  laser  driver  output  and  the  laser  input  causing  signal  reflections,  effecting 
system  performance.  Thus,  feedback  can  be  fatal  if  not  properly  taken  into 
consideration  while  designing  the  system.  So,  the  study  suggest  that  any  sys¬ 
tem  design  must  take  feedback  into  account  which  indirectly  implies  that  the 
tool  used  for  design  should  be  capable  of  including  feedback  in  system  design. 


Fig.  27.  Effect  of  increasing  feedback  on  terminal  voltage  of  laser. 
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4.6  Conclusion 


We  present  a  study  of  free  space  optic  al  system  using  fREEDA,  Our  results 
demonstrate  its  capabilities  in  understanding  the  interactions  that  can  ef¬ 
fect  the  system  performance  under  different  scenarios.  fREEDA  can  easily 
implement  any  complex  device  such  as  multi-mode  VCSEL  using  simple  set 
of  ecpiations,  whereas  same  device  modeling  in  other  circuit  simulator  would 
require  converting  these  equations  into  equivalent  circuit.  fREEDA  can  also 
capture  the  impact  of  optical  feedback  in  system  performance.  Results  show 
that  one  can  design  a  full  system  using  simple  equations  and  understand  the 
performance  related  issues  before  actually  fabricating  the  system. 
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